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Abstract. We study Hamiltonians with point interactions in spaces of vector- 
valued functions. Using some information from the theory of quantum graphs 
we describe a class of the operators which can be reduced to the direct sum 
\f} ' of several one-dimensional problems. It shown that such reduction is closely 

f— *^ ' connected with the invariance under channel permutations. Examples are 

provided by some "model" interactions, in particular, the so-called <5-, <5'-, and 
, the Kirchhoff couplings. 

<D ■ 

m : 

1. Introduction 

Quantum-mechanical Hamiltonians in coupled channels arc a natural gcncraliza- 
tion of the Schrodinger operators with zero-range interactions [1]. They provide a 
simple example of matrix-valued Fermi pseudopotentials [2] and can be described 
using the tools of the extension theory [3,4]. At the same time, such systems can be 
<^> viewed as special quantum graphs [5-8], so that one can use the general technique 

^sO ■ in order to describe all possible interactions in channels, their spectral characteris- 

tics etc. In the present paper, we develop the formalism of coupled channels using 
the quantum graph approach. 

Such an approach gives a possibility to use the self-adjoint extension theory and 
i-£h ■ the symplectic technique [9] in order to describe all possible boundary conditions. 

This can be done in many ways including the transfer matrix formalism, which 
<~| is widely used in scalar one-dimensional point interactions [10]. Using the rep- 

resentations obtained we show that a wide class of the Hamiltonians in question 
admits decoupling, i.e. by a certain unitary transformation one can reduce them 
to the direct sum of well-studied Schrodinger operators with point interactions; we 
call such Hamiltonians as well as the corresponding boundary conditions reducible. 
Such conditions can be formulated in various terms, including continuity proper- 
ties of functions from the domain of the Hamiltonian. An essential feature of the 
matrix Hamiltonians considered as quantum graphs is the presence of isometric 
parts (channels). In general, we show that the possibility of the reduction to scalar 
problems is always connected with certain invariance properties of the boundary 
conditions with respect to channel permutations. More precisely, it is proved that 
the matrix Hamiltonian is reducible iff the boundary conditions are invariant under 
the cyclic coordinate shift in a certain orthonormal basis; similar correspondence 
was found recently in connection with the inverse scattering problem on graphs [11]. 
Although such a decoupling is not a generic property, many "standard" boundary 
conditions appear to be reducible, in particular, the so-called (5, 8 P , S' , and S' s cou- 
plings as well as Kirchhoff 's boundary conditions [12] are reducible (in our opinion, 
this may illustrate the difference between the general quantum graphs and the cou- 
pled channels: the model interactions on graphs appear to be trivial in channels, 
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although one can use the same technique for their study). The reduction permits 
us to describe the spectrum of the simplest matrix Kronig-Penney Hamiltonians 
(periodically coupled channels) and to show how their parameters influence various 
spectral effects like embedded eigenvalues or the number of gaps. 

2. Parameterization of boundary conditions in coupled channels 

Consider a free particle on a graph, with the Hamiltonian acting on each edge 
as ipj i— » —ipj, where j is the edge index. Assume that the graph has the simplest 
structure, i.e. consists of n half-lines [0, oo) coupled at the origin. The boundary 
conditions take the form Aip(0) = Bip'(0), where A and B are n x n matrices 
satisfying the following two conditions [7] : 

AB*=BA*, (la) 
the block matrix (AB) has maximal rank. (lb) 

The condition l|lb(l may be rewritten in an equivalent form det(AA* + BB*) ^ 
or det(B±iA) ^ [13]. 

Now let us consider the free motion on n lines coupled at some point 5 £ i 
The Hamiltonian of the problem is the operator H — —d 2 /dx 2 acting in the space 
i 2 (K, C n ), and the coupling is described by some boundary conditions at q. For- 
mally one can consider the n coupled lines as 2n coupled half-lines, so that all 
possible boundary conditions take the form 

a(±!)=bH£A, (2) 



(3) 



J(q+)J \f'(q- 

where A and B are 2n x 2n matrices satisfying Jy. From the other point of view, 
the nature of coupled channels requests other types of parameterization [4] , namely, 
the transfer matrix formalism, 

'f'(Q+T\ _ (Cn C l2 \ (f{q- 
J(q+)J \C 21 C 22 )\f{q-) 

Below we will use mainly boundary conditions of the form J2J). Nevertheless, in 
many situations it is useful to know the connection between these two types of 
parameterization. The following proposition generalizes a construction of [3]. 

Proposition 1. The boundary conditions (j3Jl define a self-adjoint operator in 
i 2 (R, C") iff the matrices Cjk, j, k = 1,2, obey 

Ci 2 C n — C\\C\ 2 = 0, C 2 \C 22 — C 22 C 21 — 0, 

( -'ll ( -'22 ~ W2<-^21 — ^n- 

The conditions (@J are equivalent to 

^1]V21 — ^21 ( - / H — U i (^12^22 — C- 22 Ci2 — U, 
^W^ll — ^21^12 — tlin 

Proof. Substituting the equalities 

(f'(q+)\ (0 \ (f(q-)\ (0 E n \ (-f(q 
[f (<!+)) ~ \0 E n J {f(q+)J + {0 0j{ f'(qA 

'f'{Q-)\ = (0 0\ ff(q-)\ (E n 0\ f-f'(q 
f(q-)J {E n 0) {f(q+)J { 0j{ f'(q-\ 

into © we obtain 

'C12 \ (f(q-)\ _ (C u E n \ (-f'(q-) 
y C 22 -E n J \f(q+)J \C 21 J { f(q+) 



(6) 
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These boundary conditions define a self-adjoint operator iff the conditions (JTJ are 
satisfied. Eq. (|lb(l holds due to the presence of the blocks with E n , and Eq. flajl 
takes the form 

C12C11 C\ 2 C 21 \ _ / C\\C\ 2 C11C22 — 1 

K C 22 C n — 1 C 22 C 2 \) \p2lC*2 C 2 \C 22 

which is exactly These conditions means the equality 

Cn C12 \ _ / C22 — C12 
C 2 i C22/ V - ^! C11 



— ^2n, 



which is equivalent to 



-C|i C*! y VC21 C22 



#2 



and results in (0). □ 

If n = 1 (i.e. we have just one channel), the conditions are well known [10]. 
The blocks Cjk are just complex numbers, and the conditions C\ 2 C\\ = C\\C\ 2 and 
C21C22 = C22C21 mean that argCn = argCi2 = : 9i and argC 2 i = argC 2 2 =: 8 2 . 
Put C n = ae Wl , C 12 = be Wl , C 21 = ce^ 2 , and C" 22 = de i9 \ where a, b, c, d € K and 
6»i, 6*2 G [0,2tt). The third condition in Q reads as (ad - bc)e^ 02 ~ ^ = 1, which 
means that 9\ = 2 =: 0, and the boundary conditions take the form 

(/(^)) =ei9 (c d) (/(?-))' 0e[ °' 2 ^' a ' 5 ' c ' deR ' «rf-^=l. 

It is reasonable to call boundary conditions admitting the representation 
connecting. Clearly, some boundary conditions are not connecting, for example, 
the direct sum of the Dirichlet at q— and q+. Below we discuss some less obvious 
examples. 

Let us return to the parameterization J5J). The use of the values 

IV = (/(<?-), /(<?+)), r 2 / = (-/'(«-), /'(</+)) (7) 

has its origin in the theory of self-adjoint extensions of symmetric operators [9, 14, 
15]. We recall briefly some notions from the theory of abstract boundary values [16]. 
Let S be a symmetric operator in a certain Hilbert space with the domain dom S, 
S* be its adjoint with the domain dom S* . Let V be some auxiliary Hilbert space, 
and T±, T 2 be linear maps from domS 1 * to V such that 

(f,S*g)-(S*f,g) = (T 1 f,T 2 g)-(T 2 f,T 1 g), for any f,g€domS*, (8) 

and for any (vi, v 2 ) S V x V there exists / £ domS* with Tif = i>i, T 2 f = v 2 . If 
/ G domS 1 *, the values Tif and T2/ are called boundary values of /, and the triple 
(V,ri,r2) is called a boundary triple of the operator S. Boundary triple exists iff 
S has equal deficiency indices (i.e. has self-adjoint extensions), and in this case 
the dimension of V coincides with this deficiency index, see Chapter 3 in [16] for 
detailed discussion. If the space V is finite-dimensional (i.e. if the deficiency indices 
of S are finite), then all self-adjoint extensions of S are restrictions of S* on the 
elements / S dom 5* satisfying abstract boundary conditions ATif — BT 2 f, where 
A and B are matrices satisfying To obtain a one-to-one parameterization of 
the self-adjoint extensions one can normalize A and B by choosing unitary matrix 
U with 

A = l-U, B = i{l + U). (9) 

Unitary 2n x 2n matrices form a 4n 2 -dimensional real manifold, which is exactly 
the number of real parameters in the problem. 

Let q £ K be fixed. Denote by S the operator acting in L 2 (R, C n ) as ~d 2 /dx 2 
on functions from C§°(M. \ {q},C™); this operator is symmetric and has deficiency 
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indices (2n, 2n). The adjoint operator S* acts outside q in the same way on func- 
tions from VF 2,2 (R \ {<?},C™), so that the usual integration by parts in JHJ) leads 
to V — C 2 " and Y\, T 2 in the form 0, see [7]. The unitary matrix U in JHJ is 
particularly useful in approximation problems [12], and we will actively use the 
representation 10 for the boundary conditions @ . The choice of a boundary triple 
is not unique: the dimension of V is invariant, so one can always assume V = C 2n , 
and starting with given boundary operators T\, T 2 one can describe all possible 
boundary triples by means of suitable linear transformations [17]. From the point 
of view of spectral problems it may be reasonable to take as a boundary triple for 
S the set (V,fi,f 2 ) with 



ru/ 

fl 2 / 



, /(<?+) - /(<?-) / ' 2 v r22/ 



{ f(q-) + f(q+) \ 
2 

f(q-) + f(q+) 



V 



so that all possible self-adjoint boundary conditions at q take the form 



7 

(10) 



LT.f = MT 2 f, 



(11) 



where L, M are matrices satisfying the same conditions as A and B in JQl, respec- 
tively. Denote the corresponding Hamiltonian by H L ' . If the boundary conditions 
(J2J) and illlll are equivalent, then one can choose the corresponding pairs of matrices 
(A, B) and (L, M) in such a way that they satisfy 



L=h K AD 1 + BD 2 ), 
M = BDi - AD 2 , 
-E„ 



Do 



E n 

E n 
E n 



and 



A = LK 2 -hlK 1 , 
B = ~MK 2 + LK 1 , 



Ki 



E n E n 





— E n E n 



(12) 



We emphasize that due to the non-uniqueness of the parameterization this corre- 
spondence is not unique. 

For the sake of completeness we describe also the resolvents of Hamiltonians in 
coupled channels, which is useful in spectral problems. Let qi, . . . , q m , m € N, be 
points of R, qi < ■ ■ ■ < q m . Consider the operator in L 2 (M,C") acting as —d 2 /dx 2 
on functions / <E W 2 ' 2 (R \ {q : , . . . , q m }, C m ) satisfying 



/'(«.-) 
/(«.+) 



f'(Qs 
f(q s - 



/= 2 



1 ( f(q s ~) + /(?.+) 



where for each s the 2n x 2n matrices L^ s ^ and M^ s ) satisfy the same condi- 
tions as A and B before. Denote by L and M the 2mn x 2mn block matrices 
diag(L w , . . . , L im) ) and diag(M< 1 ), . . . , M (m )), respectively. The Hamiltonian de- 
scribed will be denoted by H L M . Note that the operator H° = H E ™'° is just the 
free Laplacian. The following proposition is a variant of the Krein resolvent formula 
expressed in terms of boundary conditions [18]. 
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Proposition 2. For £ € C\R+ denote by Q(C) the 2mn x 2mn matrix consisting 
of the 2n x 2n blocks Q ( '' s) (C)> 

n a.),rt e-V-^'-g-l ( ~T^ E n sign( 9! - q s )E n \ 

\-sign(<j/ - <j s ).E„ -\f<E n ) 

here and below we assume that sign = and that the square root branch is chosen 
by the condition 3V — C > f or C € (0> +°°)- ^/ swc/i ( is a regular value of H L M , 
then the matrix MQ{C) — L is invertible and for any f — (/i, ...,/„) € L 2 (R, C") 
i/ie following relation holds: 

(H L ' M -0- 1 f = (H°-()- 1 f 



~ ( H «2n( S -l)+j,2n(i-l) + fc(C)<4 '/fe)-9c S)e .7' 

s,i=l M,fc=l 
n 

- ^ a2n(s-l)+j,2n(i-l)+n+fc(C)(^i /fc)^ e i 

n 

— ^ «2n(s-l)+n+j,2n(i-l)+fe(C)(5 r ^) fk)h^ ej 



n , 
^ Q 2n(>-l)+n+j,2n(2-l)+7i+fc(C)(^ ) /fc ) ^ e j j 7 (13) 

where the numbers ctjk(C)> J, = 1, - - ■ , 2mn, are the entries of (MQ(£) — L) X A/ 7 

5 w (a; ) = _L_ e -v^i^i, 

^ W = Sign(X 2 - gs) e -v^l^l s = l,...,m, 
and 6j, j = 1, . . . , n, is the standard basis of C n . 

Using i|12|) and © one can easily express the resolvent in terms of the parameters 
A and B in J5J or C jk , j, k = 1,2, in 

3. Decoupling of the single-vertex graph 

To emphasize the specifics of the problems with coupled channels let us re- 
turn to the case of n half-lines coupled at the origin. The Hamiltonian of the 
problem is —d 2 /dx 2 acting in L 2 ((0, oo), C") = 0™ =1 L 2 (O, +oo) on functions 
/ € W 2 ' 2 ((0,oo),C n ) satisfying the boundary conditions Af{0) = Bf'(0) with 
suitable A and B from We normalize A and B by choosing them in the form 
© with suitable U € U(n); denote the corresponding Hamiltonian by Hjj- Choose 
G U(n) such that the matrix Q~ 1 UQ is diagonal. Denote by the same letter 
the associated unitary transformation of L 2 ((0, oo), C") , (0/)(x) = Qf(x). For 
x £ (0, oo) there holds (<df)"(x) = Qf"(x). This means that reduces the bound- 
ary conditions to a direct sum, 

diag(l - e w \ . . . , 1 - e* e ") 5 (0) = i diag(l + e i9 \ . . . , 1 + e t9 ")g'(0), g = 0/, 

where e lSj , j = 1, . . . , n, are the eigenvalues of U. In other words, the operator Hrj 
appears to be unitarily equivalent to the direct sum ©™ =1 -ffj, where each Hj is a 
self-adjoint operator in L 2 (0,oo) acting as —d 2 /dx 2 on functions gj € W 2,2 (0, oo) 
satisfying the boundary conditions 

(l-e"O0j(O)=i(l + e <9 Ofl$-(O), i = l,...,n. 
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Therefore, the spectral properties of Hjj are determined by the eigenvalues of U 
only, i.e. by n parameters from § ; moreover, two such n-tuples differing only by the 
order of terms are equivalent. This means, in particular, that the inverse scattering 
problem on the single vertex graph has multiple solution; note that matrix U can 
be still uniquely recovered from the scattering data [7]. 

Although the above schema gives a complete result, its applicability is rather 
restricted if the graph contains more than one vertex. One can find some gener- 
alizations for "star-shaped" graphs, i.e. if instead of the Hilbert space ©L 2 (0,oo) 
one deals with the space ©L 2 (G), where G is some graph; this models n identical 
graphs G coupled at a certain point. But even in this case the transformation O 
mentioned above is non-local and leads in general to non-local boundary conditions 
at other vertices of the partial graph G. To obtain a reasonable gain from such 
a procedure one should consider only diagonalizing transformations preserving the 
structure of G. We illustrate such a possibility by problems with coupled channels. 



4. Decoupling of channels 

Denote TL = ffi™ =1 L 2 (M) = £ 2 (K,C n ). Let Q be a uniformly discrete subset of 
K, i.e. 

inf \p - q\ = d > 0. 

On the domain domS* = Cq°(R \ Q,C n ) consider the operator S — —d 2 /dx 2 ; the 
adjoint operator S* acts in the same way on the domain domS* = W 2,2 (R\Q, C n ). 
To obtain self-adjoint operators one should introduce boundary conditions at all 
points of Q as described in Section [21 the uniform discreteness of Q guarantees that 
the operator obtained is self-adjoint [8]. Such an operator can be interpreted as the 
Hamiltonian of a free particle in n channels coupled at the points of Q. We say that 
such a Hamiltonian H is reducible iff there exists O € hi (n) such that the unitary 
transformation TL 3 f t— ► 0/ 6 TL reduces H to a direct sum of n one-dimensional 
point interaction Hamiltonians. 

All possible boundary conditions at q e Q have the form 

(1 - C/(g))r 1 / = i(l + U(q))T 2 f (Ti - tT 3 )f - U(q)(T 1 + tT 3 )f, 

ri/ = (/(<?-), /(<?+)), r 2 /= (-/'(<?-),/(<?+)), u( q )&u(2n), 141 

Denote the Hamiltonian corresponding to these boundary conditions by Hq^jj. 
Clearly, in the case of finite Q all possible Hamiltonians with point interactions 
are parameterized by 4n 2 |Q| real parameters. 
Representing U(q) in the block form, 

TT(a\ - Ul2 ^ 
[q, ~\U 21 (q) U 22 (q) 

we conclude that Hq^u is reducible if and only if the nxn blocks Ujk(q), j,k= 1,2, 
q G Q, can be diagonalized simultaneously in some orthogonal basis, i.e. if there 
exists 8 € U (n) with 

O- 1 \ /E7ii(g) U 12 (q)\ (Q 0\ (A n {q) A 12 (q) 
6"V \U 21 {q) U 22 {q)) \0 e) \A 2 i(q) A 22 {q) 

A jk(q) =diag (X jk (q,s)), 

Xjk{q, s) are the eigenvalues of Ujk(q), s = 1, . . . , n, j, k = 1,2, q £ Q. 



(15) 



The unitary transformation TL 9 / h- >• 0/ e TL reduces Hq,v to the direct sum 
of one-dimensional Hamiltonians, as due to l|15l) and to the equalities (0/)(g±) = 
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9(/(g±)) and (Qf)'(q±) = G(f(q±)) the boundary conditions (JUJ for g = 6/ 
take the form 

( 9s(q-) + i9s(q-)\ _ ( An(g, s) Xi 2 (q, s)\ ( 9s{q-~) - ig' s (q-)\ , 
\9s(<l+)- i 9s(q+)) W(q,s) \ 2 2(q > s))\g s (q+)+ig' a (q+))> s ~^---> n - 

(16) 

For the generic interaction, the condition for the reducibility of the boundary con- 
ditions can be formulated as follows: 

Proposition 3. The boundary conditions (|14Jl are reducible if and only if all the 
blocks Ujk(q), j, k = 1, 2, q G Q, are normal and commute with each other. 

This means that the reducible boundary conditions are parameterized, roughly 
speaking, by n\Q\ unitary 2x2 matrices \\jk(q, s)) . fe=1 2 , q £ Q, s = 1, . . . , n (up 
to permutations), and a unitary n x n matrix which diagonalizes the boundary 
conditions. This means that the most general reducible boundary conditions involve 
n\Q\ dimuW(2) + dim^ U{n) = 4n\Q\ + n 2 real parameters. 

An analogue of Proposition |31 can be given in terms of the transfer matrix 0. 

Proposition 4. The Hamiltonian Hqjj given by the boundary conditions © is 
reducible iff the blocks Cjk, J, k = 1,2, are normal and commute for all q £ Q. In 
particular, if Q consists of a single point q and Cjk are self-adjoint, then Hq jj is 
reducible. 

Proof. The first part is obvious. Assuming that that the blocks Cjk, j,k = 1,2, 
are self-adjoint (like it was done in [4]), one concludes from Q and 10 that they 
all commute with each other and, therefore, can be diagonalized simultaneously. 
But this means that the corresponding Hamiltonian is reducible. □ 

It is useful also to have "quantitative" reducibility criteria in terms of bound- 
ary conditions. The corresponding matrix U may be difficult to find, but the 
reducibility can be found by other means. To illustrate this, we consider the 
Hamiltonian H given by its quadratic form Q{f,f) = (/',/') + (f(q),Af(q)}, 
domQ = W 1,2 (R,C n ), where A is a n x n self-adjoint matrix. This Hamiltonian 
may be viewed as the so-called matrix (S-potential and corresponds to the boundary 
conditions /(<?-) = f(q+) =: f(q), f'(q+)-f(q-) - Af(q), cf. [8]. Clearly, an or- 
thogonal transformation which diagonalizes A will reduce the boundary conditions 
to a direct sum; the Hamiltonian H is unitarily equivalent to the operator —d 2 /dx 2 
with the boundary conditions gj(q-) = gj(q+) -■ gj(q), 9j(q+)-9j(q-) = ctjgj(q), 
j = 1, . . . , n, respectively, where ctj are the eigenvalues of A. In other words, H is 
isomorph to the direct sum of the usual one-dimensional ^-perturbations. Let us 
try to generalize this example. 

Propositions. LetQ consist of a single point q . If there exist a, (3 £ M., \a\-\-\f3\ > 
0, and c, c' £ { — 1, 1} such that 

a(f'(q+) + c'f'(q-)) = 0(f(q+) + cf(q-)) (17) 

for all f £ doTaHjj , then Hjj = Hqjj is reducible. 

Proof. Consider first the case c = c' = — 1. Assume first a = 0. Put D = {/ £ 
dom^* : f(q—) = f(q+) ='■ /(?)}• Clearly, dom Hu C D, and for arbitrary 
f,g £ D there holds 

0V, r^) - (r 2 /,r l5 ) = -</(?-), g' (<?-)) + (/(?+),</(<?+) 

+ </'(?-), <?(<?-)> - (/'(<?+), 9(q+)) = Kf,r' 2 g) - (n/,r; ff ), 

where T'J = f(q), T' 2 f = f(q+) - f(q~). Denote by S the restriction of S* to 
the set dom 5*o = {/ £ D : T'J = T' 2 f = 0} = {/ £ dom S* : /(<?-) = f(q+) = 
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0, f'(q+) = /'(?— )}< Clearly, this is a symmetric operator, and the set D is the 
domain of its adjoint Sq. Therefore, (C n ,r' 1 , L 2 ) is a boundary triple for this new 
operator So- As Hjj is a self-adjoint extension of So, there exists V G U(n) so 
that Hjj is determined by the boundary conditions (T[ — iT' 2 )f — V(T[ + iT' 2 )f, 
f G dom^Q. Let O be a unitary transformation which diagonalizes V. Clearly, 
induces a unitary transformation of H., and the components of the function g = O/, 
/ G dom Hjj, satisfy 

#(?-) = aM+) =■ g(q), (i - ^ = i(i + ^) • (g'j( q +) 9j ( q -)), (is) 

e l6,j are eigenvalues of V, j = 1, . . . , n. (19) 

Therefore, Hjj is reducible. 

Consider now the case a ^ 0. Put 7 = /3/<x We use the boundary triple (|l(Jf> . 
Denote by D the set {/ G dom S* : f u / - jt 12 f = 0}. The condition l(T7|) 
means the inclusion domii)y C -D. Let /, g G -D, then (Fx/, L 2 g) — (L 2 /, Lig) = 
(Fi/,r 2 . 9 ) - (L 2 /,F' l3 ) with T'J = f 12 /, T' 2 f = 7 f 21 / + f 22 /. Denote by 5 the 
restriction of S* to the set dom So = {/ € D : T^f = T' 2 f = 0}; this is a symmetric 
operator, D — domSg, and (CT'Tj) is a boundary triple for So- As Hjj is a 
self-adjoint extension of So, there exists V £U(n) such that Hjj is determined by 
the boundary conditions (Li - iT' 2 )f = V(T[ + iT' 2 )f, f G D. Let 6 G U{n) such 
that <d~ 1 V@ is diagonal. Noting that 9 commutes with all the operators Tjk, Tj, 
j, k = 1, 2, we reduce the boundary conditions to a direct sum for g = Of. 

Now let d = — 1, c = 1. Denote by D the set of functions / G domS 1 * satisfying 
HI 70 and use again the boundary triple (JTDJ, then for any /, g G -D there holds 
(Li/,L 25 ) - (L 2 /,L l3 ) = (Li/,L 2 . 9 ) - (L 2 /,L' l3 ) with T'J = L 12 /, L 2 / = L 22 /. 
Denote by So the symmetric operator which is the restriction of S* to the domain 
dom S = {/ G D : T' 1 f = T' 2 f = 0}, then D = domSp. Taking into account the 
fact that Hu is a self-adjoint extension of So we proceed with the proof as in the 
previous case. The rest combinations of c and d can be considered in the same 
way. □ 

To formulate an important corollary we recall that a function / : M — ► M is called 
anticontinuous at q if there exist the limits f{q±) and f(q+) + f(q—) = 0. 

Corollary 6. Let the set Q consist of a single point q. If one of the following 
conditions is satisfied: 

• all functions from AoyuHq^jj are continuous, 

• all functions from dom Hq^u are anticontinuous, 

• derivatives of all functions from dom Hq^jj are continuous, 

• derivatives of all functions from dom Hqjj are anticontinuous, 

then Hqjj is reducible. 

We emphasize again that the last proposition and the corollary apply to chan- 
nels coupled at one point only. Of course, this works also for channels which are 
identically coupled at several points. 

5. Permutation-invariant boundary conditions 

Let us return to the general Hamiltonian Hqjj with boundary conditions at 
point of a discrete set Q (see the beginning of the previous section) . The aim of 
this section is to discuss a correspondence between the reducibility and invariance 
under channel permutations. 

The most general version of this correspondence can be formulated as follows: 



REDUCIBLE BOUNDARY CONDITIONS IN COUPLED CHANNELS 







Proposition 7. A matrix-valued point interaction Hamiltonian Hqjj with a point 
interaction supported by a uniformly discrete set Q is reducible if and only if there 
exists a unitary n x n matrix with non- degenerate eigenvalues such that the 
boundary conditions at all points of Q are invariant under the transformation f i— ► 



Proof. At each point q £ Q there exists U = U(q) <E U(2n) such that all the 
functions from the domain of Hqjj are characterized by the condition JHJ. 

If the Hamiltonian is reducible, all the blocks Ujk(q), j,k = 1,2, q G Q, are 
diagonal in some orthogonal basis and. therefore, commute with any matrix which 
is diagonal in this basis. Taking an arbitrary diagonal unitary matrix with non- 
degenerate eigenvalues we show that the condition formulated in the proposition is 
necessary. Let us show that this condition is also sufficient. 

The invariance of boundary conditions under means that all the blocks Ujk(q), 
j,k = 1, . . . , 7i, commute with 0. This means that the invariant subspaces of are 
such for all blocks at all points of Q. As these subspaces are one-dimensional and 
orthogonal to each other, all the blocks are diagonal in the eigenbasis of 0. □ 

This proposition shows that the reducibility is an effect which is closely connected 
with non- uniqueness in the inverse scattering or spectral problems on graphs [11]: If 
the Hamiltonian is invariant under a unitary transformation with certain properties, 
then there exists another graph (in our case, the union of real lines with marked 
points) having the same spectrum. 

An important example is provided by Hamiltonians which are invariant under 
certain channel permutations. 

Corollary 8. Let a be a permutation of order n (i.e. a n = id and o~ k ^ id for all 
k £ {1, . . . , n — 1}) and Hqjj be invariant under the transformation fj t— > f a ij), 
j = 1, . . . , n, then Hq jj is reducible. 

Proof. Indeed, the minimal polynom of the transformation is A™ — 1, which means 
that all eigenvalues are simple. □ 

Actually, this situation is in some sense generic, as the following proposition 
shows. 

Proposition 9. The Hamiltonian Hqjj is reducible iff there exists an orthonormal 
basis (hi, . . . , h n ) in C n , so that all boundary conditions are invariant under the 
transformation hj i— > h^_ 1 - )modn . 

Proof. The minimal polynom of the transformation described is again A™ — 1, 
which means that all eigenvalues arc simple. This shows that the existence of such 
transformation is sufficient for the boundary conditions to be reducible. Let us 
show that this condition is also necessary. 

Let Hqjj be reducible, then there exists an orthonormal basis G — (gj,j = 
1, . . . , n) in C™ in which all blocks Ujk(q) are diagonal. In this basis, define a linear 
transformation 3 by its matrix diag(Ai, . . . , A„), Aj = exp(27rij/n), j = 1, . . . ,n. 
Clearly, the blocks Ujk{q) commute with 2 (as all these matrices are diagonal). 
From the other side, in the basis 



e/. 
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the transformation S has the matrix 

/ 1 
1 




V 1 



0\ 


1 



which is exactly the matrix of the transformation 



(hj,f)hj i ^ ^{hjj)^-! 

Therefore, the boundary conditions at all points are invariant under the cyclic shift 
of coordinates with respect to the basis (hi, ... , h n ). □ 



Of certain interest are Hamiltonians (and the corresponding boundary condi- 
tions) which are invariant under all channel permutations. Clearly, this means that 
the blocks Ujk{q) are of the form Ujk(q) = a,j k (q)E n + bj k (q)J n , where J„ is the 
nx n matrix whose all entries are equal to 1, and the complex numbers aj k {q) and 
bjk(q) obey the condition (a jk (q)) j fc=12 , (a jk (q) +nb jk (q)) j k=12 & U{2). (Clearly, 
the spectrum of J n consists of a simple eigenvalue n and a (n — l)-fold degener- 
ate eigenvalue 0.) This class includes the frequently used <5, S' s , 8 pi and 5' cou- 
plings, which we consider in greater detail (some different notation is used, see [2]). 
For more detailed discussion of the origin of these coupling types we refer to the 
works [12,19] and references therein. The corresponding boundary conditions for a 
function / e W 2 ' 2 (R \ {q}, C") are as follows: 



5(q,a) 



5' s (q,P) 



S p (q,a) 



m-) = fk(q+) 

n 

£ (fM) -fte-)) 

-f'M-) = /*(«+) 

n 

E(/i(«+) + /iH 

±f j (q±)Tf' k (q±) 

n 

£(/,(<h +/.(<?+)) 



j,k = 1, . . . , n, 



(20a) 



(20b) 



fj(q±) - f k (q±) 



f(q), j,k = i, 

af(q), 

/'(«), 
Pf(q), 

(fj(<l±) - fk(q±)), j,k = i,...,n, 
o. 

(20c) 

±/j(«±)T/fc(«±)), j,fc = l,...,n, 



a 
2^ 



5' (<?,/?) 



A, 

2n 



E(/i(9+)-/^-)) = ' 



i=i 



(20d) 



where a and /3 are real parameters. The 5(q, 0)-coupling corresponds to the so-called 
Kirchhoff boundary conditions at q; they appear, for example, if one considers the 
coupled channels as a limit of shrinking manifolds [20] . For the sake of brevity we 
denote the introducing of boundary conditions as a formal sum, for example, under 
the operator 

H = -^+6' s (q u p)+6(q 2 ,a) (21) 
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we mean the operator which acts as / i— > — /" on functions / G W 2,2 (M\{qi, q 2 }, C n ) 
satisfying the boundary condition l|20b(l for q = qi and Ij20a|) for q = q 2 . In onc- 
dimcnsional case we use a more traditional way of writing, for example, 

+ /38' s (x - qi ) + aS(a - q 2 ) (22) 



H = — 



dx 2 



will denote the same operator as in (|21() assuming that n = 1. In fact, one can 
consider the expression ill' 21) as a self-adjoint operator if one uses the theory of 
distributions with discontinuous test functions [21], see also [2]. 

Proposition 10. Let Q, Q' s , Q p , Q' be non-intersecting discrete subsets of K, 
and their union P := Q U Q' s U Q p U Q' be uniformly discrete. Denote by H the 
self-adjoint operator in L 2 (R, C n ), n > 1, of the form 

d 2 



q£Q 



qeQp 



where a q , (3 q are real parameters. Then H is unitarily equivalent to the direct sum 
©fe = i-fffc, where H}~ are self-adjoint operators in L 2 (R), namely, 



H i = -h+H - E - *.(*-«)+ E - ^(»-?)+ E -n^-q), 

dx 2 - ' n ' n n n 

q£Q qEQ' s qEQp qEQ' 



i.e. the operator —d 2 /dx 2 acting on functions f G T / K 2,2 (R \ P) satisfying 



(23) 



/(<?-) = /(?+) =: /(?), 

/'(«-) + /'(?+) = o, 
f(q-)+ /(?+) = 0, 

/'(?-)=/'(«+) =:/'(?), 



/'(?+)-/'(?-) 
/(<?-)+/(<?+) 



/'(?+) -/'(?-) = ^ /'(«+), 

/(«+)- /(«-) = ^/'(g), 



?6Q' 



and the operators H 2l ■ ■ ■ , H n are equal to each other and act as g(x) t— > —g"(x), 
x ^ P, on functions g G W 2,2 (R \ P) satisfying the following boundary conditions: 

9{q-) = g(q+) = 0, gGQ, g'{q-) = g'(q+) = 0, q £ Q' s , 

ot q g{q~) + 2ng'(q-) = a q g(q+) - 2ng'(q+) =0, G Q p , (24) 
2n$(g-) + /3 gff '(g-) = 2ng(q+) - f3 q g'(q+) = 0, q G Q'. 

Proof. We recall that the uniform discreteness of P guarantees the self-adjointness 
of if [8]. 

As it was shown in [12], the boundary conditions 1)20(1 can be written as 1)14(1 
with U(q) = a q E 2n + b q J 2ni where 



-1 
1 

In — ia q 

2n + ia q 
2n + iJ3g 

2n — if3 q 



for 5(q, a q ), 
for S' s (q,/3 q ), 

for 5 p (q,a q ), 
far 



2n + ia 



2n — i[3q 
2 

2n + ia„ 



for 5(q, Oig), 

for S'MPi), 
for S p {q, a g ), 



(25) 



„ 2n — ip c 



J 2n for S'(q,/3 q ). 



The n x n blocks of U are of a rather simple form, namely, U\\{q) = U 22 (q) = 
a q E n + b q J n , U\ 2 {q) = U 2 ±(q) = b q J n . Let S be a linear transformation which 



12 



KONSTANTIN PANKRASHKIN 



diagonalizes J n , then at each point q € P the components of the functions g := S/, 
/ € domif, satisfy 

U(9+)-^fc(9+V ^UM+WJ' (26) 

with 

«w-(-^ at to- - 

which is exactly (j2HJ) and J23}. □ 

All the boundary conditions (|24J) are obviously non-connecting; this means that 
none of the couplings (|20|) admits the representation © . This is connected with the 
fact that these couplings are actually invariant also under half-channel permutation. 



6. Periodically coupled channels 

Let us illustrate the separability effects by periodic problems with point interac- 
tions. Periodically coupled channels provide simple examples of periodic quantum 
graphs, so that the general powerful technique for their analysis is available [22,23]. 
The previous discussion gives a possibility to describe the spectrum of some peri- 
odic Hamiltonians by other means: one can easily reduce the spectral problem for 
periodically coupled channels to the spectral problem for periodic scalar Hamilto- 
nians with point interactions, i.e. to the well-studied generalized Kronig-Pcnney 
models [24]. We restrict ourselves by considering some examples. 

Example 11 (Permutation-invariant delta-potential). In L 2 (R, C 2 ) consider the 
periodic delta-potential invariant under channel permutation; this corresponds to 
the boundary conditions 

f(q-) = f(q+)=-f(q), f(q+)-f'(q-) = (aE 2 +{3J 2 )f(q), a,0&R, q e ttZ. 

Elementary considerations show that this Hamiltonian H is unitarily equivalent to 
the direct sum Hi © H 2 , 

d 2 d 2 
Hi = -— + a^2S(x - nn), H 2 = + (a + 2/3) 2J S(x - nn), 

so that the spectrum of H is the union of the spectra of Hi and H 2 . If both a 
and a + 2/3 have the same sign, then the spectrum of H has an infinite number 
of gaps. For example, for a, a + 2(3 > the spectra of Hi and H 2 consist of the 
bands (a m ,m 2 ) and (b m ,m 2 ), m = 1,2,..., respectively, where a m ,b m > (m — 
l) 2 , see Theorem III. 2. 2. 3 in [1]. The spectrum of H consists then of the bands 
(min(a m ,6 m ),TO 2 ), m = 1,2, .... 

Let us show that H has only a finite number of gaps if a(a + 2/3) < 0. To be 
definite, assume that a > and a + 2/3 < (the second case can be considered in 
the same way). The spectrum of Hi consists of the bands (a m ,m 2 ), m= 1,2,..., 
where a m = m 2 — 2m — 2a/Tc — l + 0(l/m), m — > 00, and the spectrum of H 2 consists 
of the bands (A m , B m ), m = 1,2, ... , where A\ < Bi < 0, B m > A m = (m — l) 2 , 
to = 2,3, . . ., B m = (m - l) 2 + 2m + 2(a + 20) - 1 + 0(l/m), to -> 00, see 
Theorem III. 2. 2. 3 in [1]. Obviously, for large m there holds B m > a m , which means 
that the gaps are overlapped by the large bands. 

Example 12 (Periodic 5-coupling) . For any real a the operator 
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acting in L 2 (R, C") is unitarily equivalent to the direct sum H a ®(®^ = ^HD), where 

lei 

and Hd is the Laplace operator in I? (R) acting on functions satisfying the Dirichlct 
boundary conditions at the points irl, I e Z. Therefore, the spectrum of H consists 
of the spectrum of H a and of the infinitely degenerate eigenvalues m 2 , to e N. 

For a 0, the spectrum of H a consists of values k 2 satisfying the Kronig-Pcnncy 
equation | cos7rfc + a/(2nk) sin7r/c| < 1, 3fc > 0, and the band edges are given by 
the values k 2 with cos7r/j + a/(2nk) sin-7r/c = ±1, see [1, Theorem III. 2. 3.1]. In 
particular, the Dirichlct eigenvalues are situated on the band edges. 

Example 13 (Periodic Kirchhoff coupling). Let us emphasize a particular case 
of the previous example. If a = (Kirchhoff couplings), then H a is just the free 
Laplacian. The spectrum of the initial operator H, i.e. of channels periodically 
coupled by the Kirchhoff boundary conditions, consists of the semiaxis [0, +oo) and 
embedded Dirichlct eigenvalues to 2 , m = 1,2, ... . 

The existence of eigenvalues in the spectrum of a periodic problem on the graph 
in our toy situation is connected closely with the reducibility of the boundary 
conditions. Nevertheless, such effects appear in much more general structures [25]. 
It is known that a periodic graph can have eigenvalues only in the case of compactly 
supported solutions [22]. The existence of such solutions is possible only in the case 
of the so-called analytically disjoint couplings, which can produce even stronger 
spectral effects [26]. 
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